Unusual diffuse X-ray source in the Galactic center region 
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ABSTRACT 

We report the ASCA and Chandra discovery of a diffuse X-ray source 
in the Galactic center region. The X-ray spectrum is fitted with a non- 
equilibrium ionization (NEI) plasma model of about 6-keV temperature. 
The model requires higher than solar metal abundances, a young plasma 
age of ~ 100 years and a large ATh value of about 10 23 cm -2 . The TVh 
value constrains the source position to be in the Galactic center region 
at about 8.5 kpc distance. The high resolution X-ray image with the 
Chandra ACIS shows a ring of 10" radius which corresponds to 0.4 pc 
at the Galactic center, and a tail-like structure. Although the morphol- 
ogy is peculiar, the other X-ray features are likely to be a very young 
supernova remnant, possibly in a free expansion phase. 

Subject headings: Galaxy: center — ISM: individual (GO. 570— 0.018/CXO 
J174702.6-282733 ) - - Supernova remnants - - X-rays: ISM 
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1. Introduction 

The Galactic center (GC) region within 
~100 pc is very rich in various X-ray sources, 
such as X-ray binaries, young stellar clusters 
and supernova remnants (SNRs). A big mys- 
tery is a diffuse hot (~ 10 keV) plasma of 
1° x 2° elliptical shape with the total energy 
of about 10 54 ergs (Koyama et al. 1989, 1996; 
Yamauchi et al. 1990). The dynamical age 
of the hot plasma is estimated to be about 
10 5 years. One possible origin of the dif- 
fuse plasma is multiple supernova explosions 
(10 3 supernovae) in the past 10 5 years (Ya- 
mauchi et al. 1990), which predicts that many 
young SNRs should be discovered with a deep 
exposure above the 2 keV X-ray band. With 
the ASCA Galactic center survey project, 
we found an extended X-ray source on the 
Galactic plane near the giant molecular cloud 
Sgr B2, which exhibits unusually strong iron 
lines (Sakano et al. 1998). The Chandra Cy- 
cle 1 observation confirmed the ASCA re- 
sults and revealed a ring and tail structures 
in the iron K a line band, which we named 
G0.570-0.018/CXO J174702.6-282733 (here 
and after, GO. 570— 0.018). This paper reports 
the ASCA and Chandra results of GO.570-0.018 
and discuss the nature of this peculiar source. 

2. Observations and Data Reductions 

ASCA (Tanaka et al. 1994) observed the 
Sgr B2 region near the Galactic center on Oc- 
tober 1 in 1993 (PV phase) and September 
22-24 in 1994 (A02). ASCA was equipped 
with two Gas Imaging Spectrometers (GIS2 
and GIS3) and two Solid-state Imaging Spec- 
trometers (SIS0 and SIS1) on the focal planes 
of four X-ray telescopes (XRT) (Selemitsos 
et al. 1995; Ohashi et al. 1996; Burke 
et al. 1991). In both the observations, 



GO. 570— 0.018 was placed on the center of the 
SIS-field, hence was intercepted by the center 
gap of each CCD chip. We therefore used the 
GIS (GIS2 and GIS3) data only. For GISs, 
nominal bit-assignment of PH-mode wes used 
in high and medium bit rate. After the stan- 
dard data screening, we obtained the effective 
exposure time of PV and A02 of 17.5 ksec and 
80.0 ksec, respectively. 

GO.570-0.018 was in the ACIS-I field of 
view when Chandra observed the Sgr B2 re- 
gion on 29-30 March 2000. The satellite and 
instruments are described by Weisskopf et al. 
(1996) and Garmire et al. (2001), respec- 
tively. Data reduction is made with the same 
method of the Sgr B2 analysis (Murakami, 
Koyama, & Maeda 2001). GO.570-0.018 is 
located at the top of chip 13, and thus is 
heavily affected by the Charge Transfer In- 
efficiency (CTI). To minimize the effect of 
the CTI degradation, we use the software 
developed by Townsley et al. (2000). We 
correct the energy gain using the instrumen- 
tal emission line from Ni (K a =7.5keV) and 
Au atoms (L Q =9.7keV) with an accuracy of 
< 0.5% (90% confidence). The response ma- 
trix is made using the nearly contempora- 
neous observation of reference lines from an 
on-board calibration source (OBSID=62097), 
which are analyzed with the same CTI cor- 
rection. The effective area of the telescope 
mirrors and the detection efficiency of ACIS 
are calculated with the mkarf program in the 
Chandra Interactive Analysis of Observations 
Software (CIAO, Version 2.1). After screen- 
ing the data, we obtain the effective exposure 
time of 100 ksec. 

3. Results and Analysis 
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3.1. X-ray Morphology 

In the ASCA GIS image pointed at the 
Sgr B2 cloud, we find a local excess in the west 
of Sgr B2. This excess is clearly seen if we 
limit the X-ray energy band of 6.0—7.0 keV, 
which include the iron K-shell transition line. 
Figure 1 shows the 6.0—7.0 keV band image. 
The brightest source located at the north- 
east is the giant molecular cloud Sgr B2 (Mu- 
rakami et al. 2000). The other excess at 
the center of the image is a newly discovered 
X-ray source. Due to the limited statistics 
and spatial resolution, error of the peak posi- 
tion of the diffuse structure is rather large of 
0?55 < I < 0?57, -0?02 < b < 0?00. 

In order to investigate the accurate po- 
sition and fine spatial structure of this dif- 
fuse source, we inspect the Chandra image in 
the 4.0—7.0 keV band. As is shown in Fig- 
ure 2, we see a complex X-ray structure with 
a ring of about 10" radius and an east-to- 
west tail from the ring. The center coordinate 
of the ring structure is R.A.= 17 h 47 m 02!6, 
Dec= -28°27'33" (epoch 2000), which corre- 
sponds to (/, b) = (0?570, -0?018), the same 
position of the ASCA diffuse source within 
the error. We hence designate this source as 
G0.570-0.018/CXO J174702.6-282733 (here- 
after GO.570-0.018). To investigate the sta- 
tistical significance of the shell-like structure, 
we make a radial profile with the center at 
(17 h 47 m 02 s .6, -28°27 / 33' / ) J 2ooo- Figure 3 shows 
the radial profile in the 4.0—7.0 keV band. 
The excess around 10" is statistically signifi- 
cant, hence the ring-structure is real. On the 
other hand, a faint point-like structure at the 
shell center seen in Figure 2 is not significant 
because it contains only 3 photons. Unlike 
Sgr B2, no radio nor any other wave band 
counterpart is found for GO. 570— 0.018 (see 
Oka et al. 1998a, 1998b; Tsuboi, Handa, & 



Ukita 1999; Lis, k Carlstrom 1994; Mehlinger 
et al. 1992, 1993; Seiradakis et al. 1989). 

3.2. X-ray Spectra 

The ASCA X-ray spectrum is made using 
the data in a circle of 3-arcmin radius and 
subtracting the background data in a 3' < r < 
4' ring as is shown in Figure 1. The spec- 
trum exhibits prominent line at 6.0—7.0 keV 
and large absorption at low energy band. We 
then fit with a phenomenological model of 
a thermal bremsstrahlung plus a Gaussian 
line. The line energy is determined to be 
6.60lo;o9 keV with the equivalent width of 
3.7lf;2 keV (here and after, errors are 90% 
confidence unless otherwise noted), indicating 
a K a line from iron atoms with lower ioniza- 
tion states than helium-like (6.7 keV). The 
temperature is constrained to be higher than 
2 keV. Therefore the equilibrium ionization 
state of iron should be helium or hydrogen 
like with the line energy of 6.7—6.9 keV. We 
thus conclude that the plasma is still in an 
ionizing phase or in a non-equilibrium ioniza- 
tion (NEI). 

The Chandra X-ray spectrum is made from 
the same area as with ASCA, while the back- 
ground spectrum is taken from a source-free 
region with the same CTI effect as the source 
region. The background-subtracted flux is 
20400 counts and the spectrum shows a clear 
iron line as is already found with ASCA. In 
the diffuse source region, we also find many 
faint point-like sources. We hence execute 
the CIAO 'wavdetectf software of a wavelet 
method (Freeman et al. 2000) and resolve 18 
point sources. The integrated flux (2.0—10.0 
keV) of all the point sources is 410 counts, 
which is only ~ 2% of the diffuse X-rays. 
Therefore a contamination of the point-sources 
can be ignored practically. The Chandra spec- 
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trum is well fitted with the same phenomeno- 
logical model as used for the ASCA spec- 
trum, a bremsstrahlung plus a Gaussian line. 
The plasma temperature is constrained to be 
> 3.4 keV, while the line energy is 6.501q;oI 
keV with the equivalent width of 4.11J; 4 , keV. 
Thus the best-fit parameters are consistent 
with, and are more accurate than the ASCA 
results. 

Since both the ASCA and Chandra spec- 
tra are fitted with the same model and the 
line center energy is consistent with an NEI 
plasma, we simultaneously fit the two spec- 
tra with a more realistic NEI plasma model 
(Borkowski et al. 2001). In addition to a con- 
ventional thin thermal plasma model, an NEI 
model includes another parameter, which is 
called the ionization parameter r—nt, where 
n and t are plasma density and elapsed time 
after the plasma is heated-up. As the param- 
eter r becomes larger than 10 12 ~ 10 13 cm -3 
s, an NEI plasma approaches to a collisional 
ionization equilibrium (CIE) plasma. In ad- 
dition to the strong iron line, the Chandra 
spectrum includes emission lines at the en- 
ergies of K-shell transition of argon, calcium 
atoms. We therefore vary the abundances of 
all the elements collectively fixing the rela- 
tive ratio to be solar (Anders, & Grevesse 
1989). Free parameters are the electron tem- 
perature (kT), the abundance (Z), the ion- 
ization parameter (nt), and interstellar ab- 
sorption (Ah). For the other free parame- 
ter, we allow one normalization factor to be 
fitted both the ASCA and Chandra fluxes si- 
multaneously. The NEI fit is acceptable with 
the best-fit spectra and parameters shown in 
Figure 4 and Table 1. The data excess of 
Chandra and deficit of ASCA above and be- 
low the iron line are attributable to the gain 
uncertainty of these instruments. However 



this small systematic difference is not a se- 
rious problem in the present analysis. The 
ionization parameter of r =1.7xl0 10 cm -3 s 
is very small, while the abundance is signifi- 
cantly larger than the solar value. 

4. Discussion 

ASCA discovered a diffuse X-ray source 
(GO. 570— 0.018) with a strong iron line near 
the Sgr B2 cloud, then Chandra found a 
ring-like structure of 10"-radius. The hy- 
drogen column density is determined to be 
A H = (13.9l||) x 10 22 H cm- 2 . Since the 
A H value is nearly equal to those of the Galac- 
tic center sources of the same Galactic latti- 
tude (Sakano 2000), GO.570-0.018 would be 
located near at the Galactic center region. 
We hence assume the source distance to be 
8.5 kpc. The X-ray luminosity is then esti- 
mated to be ~ 10 34 ergs and the size of the 
ring radius is 0.4 pc. The spectra are well fit- 
ted with an NEI model with plasma temper- 
ature (kT) of 6.1 keV, ionization parameter 
(r) of 1.7x 10 10 cur 3 s, and metal abundances 
(Z) of 4.5 solar. Therefore, together with the 
ring-like morphology, GO. 570— 0.018 is likely 
a young SNR, either in an adiabatic or a 
free expansion phase. We first apply a Sedov 
self-similar model assuming in an adiabatic 
phase. This model gives the emission mea- 
sure (E.M. = n 2 V), radius (R) and electron 
temperature (kT) as follows (Ostriker et al. 
1988); 

n 2 V = (An a f x 4nR 2 x (§) 
(cm -3 ), 

R = 5.0x( T? ^^) 1 /5(^^)-l/5( *a_)2/5 
^ 10 01 ergs ' v 1 cm '' ' v 10 J years ' 

(pc), and 

tvl — 1.0A^ 1()51 erg J y lcm -'A) Uo 3 years J 

(keV), 
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where t s , n a and E are the age of the 
plasma (in unit of year), ambient density (in 
units of cm -3 ) and explosion energy of the su- 
pernova (in units of ergs), respectively. Using 
the observed values of R = 0.4 pc, kT = 6.1 
keVand-E.M. = 8.1 xlO 56 cm' 3 , we obtain t s , 
n a and E to be 70 years, 5.1 cm -3 and 3.5 x 
10 48 ergs, respectively. The swept-up mass is 
then calculated to be n a x 4/37ri? 3 ~ O.O3M , 
which is extremely smaller than that of the 
SN ejecta of a few M . The estimated explo- 
sion energy of GO. 570— 0.018 is also extremely 
small compared with a usual SNR. These in- 
dicate that only a tiny fraction of the explo- 
sion energy (~ 10 51 ergs) and ejected mass (a 
few M ) have been converted to the thermal 
plasma. These facts strongly indicate that 
GO. 570— 0.018 is not in an adiabatic phase but 
still in a free expansion phase. 

The X-ray ring structure of GO. 5 70— 0.018 is 
similar to that of SN 1987A (Burrows et al. 
2000), where a strong stellar wind from a mas- 
sive progenitor might produce a gas ring of 
sub-pc radius, and was heated by the collision 
of supernova ejecta. From the radial profile in 
Figure 3, the thickness of the ring (FWHM) is 
about 6", or 0.2 pc, hence the plasma volume 
is 1.2 xlO 55 cm 3 . Since the X-ray flux from 
the ring is 67% of that of the whole area of 
GO. 570-0. 018, the emission measure (E.M.) 
from the ring is 5.5 xlO 56 cm -3 . Then us- 
ing the best-fit ionization parameter (nt) of 
1.7xl0 10 cm -3 s, the plasma density (n) and 
ionization age (t) of the X-ray ring are esti- 
mated to be 6.7 cm' 3 and 80 year, respec- 
tively. In Table 2, we compare the physi- 
cal parameters of the ring of GO. 570— 0.018 
with those of SN 1987A. The ring size and 
age are about 3—6 times larger than those of 
SN 1987A. Therefore GO.570-0.018 would be 



a "future SN 1987A" after a free expansion 
with the expected speed of v — R/t ~ 4900 
km s _1 . 

One problem of this scenario for GO. 570— 0.018 
is an additional X-ray structure, the east-to- 
west X-ray tail. We suspect that the structure 
would be made by some instability of dense 
stellar wind from progenitor, and was heated- 
up by the collision of the SN ejecta. In this 
case, the front of free expanding ejecta would 
be about two times larger than the ring ra- 
dius, hence the age of SNR is doubled. This 
age is still in a very young phase, thus the 
above discussion based on a free expansion 
phase is essentially not changed. 

In the young SNR scenario, we expect non- 
thermal radio emissions associated with the 
thermal X-rays. However we see no significant 
radio continuum flux from GO. 570— 0.018. In 
general, the X-ray flux is proportional to the 
square of electron density (n 2 ), while the ra- 
dio flux is proportional to the square of mag- 
netic field, B 2 . Assuming that magnetic field 
is proportional to the ambient gas density 
(the magnetic field is frozen to the ambient 
gas), the flux ratio between X-ray to radio 
is roughly constant among shell-like SNRs. 
Referring the radio and X-ray data of typ- 
ical young shell-like SNRs, Cas A, Tycho 
and Kepler, we estimated the radio flux of 
GO.570-0.018 to be 0.1-1 Jy (except Cas A) 
at 1 GHz, using the X-ray flux of ~10~ 12 ergs 
s _1 . SN 1987A has almost the same X-ray flux 
as that of GO. 570— 0.018, and the radio flux 
density is ~40 mJy at 4.7 GHz (Burrows et 
al. 2000). These radio flux density would be 
below the current detection limit of the radio 
SNRs located near the GC, where the radio 
background is extremely high (~0. 1—0.5 Jy). 
Thus we encourage more deep and fine spatial 
resolution search for radio emissions near at 
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GO.570-0.018. 
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Fig. I. — ASCA GIS image in the 6.0—7.0 keV (iron line) band. The dotted lines are the Galactic 
coordinate with the grid spacing of 5'. The source at the lower-right is GO. 570— 0.018 and the 
spectrum is taken from the solid circle of a 3'-radius. Background spectrum is taken from the 
annulus defined by the solid and dashed (radius = 4') circles. The brightest source in the upper- 
left is the Sgr B2 cloud. 
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Fig. 2.— Chandra ACIS-I image of GO.570-0.018 in the 4.0-7.0 keV band. The dotted lines are 
the Galactic coordinate with the grid spacing of 20". GO. 570— 0.018 is appeared to be a 10" radius 
ring with an east-west tail. This image is smoothed with a Gaussian of a — 2 pixels. Before the 
smoothing, the brightest pixel has 3 photons. The brightness increases with a linear scale from 
0.1 (black pixel) to 1.3 photons (white pixel). 
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Fig. 3. — The radial profile of the ring structure of GO. 570— 0.018. The ring center is (17 h 47 m 02!6, 
— 28°27 / 33")j2ooo- 
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Fig. 4. — The X-ray spectra of GO. 570— 0.018 obtained from ASCA (gray) and Chandra (black). 
The data points are given by crosses while the solid lines are the best-fit NEI plasma model. The 
data residuals are shown in the lower panel. 
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Table 1 

Best-fit parameters of NEI model obtained from the combined fit of the 

ASCA and Chandra spectra 



Wh 1 
[10 22 Hcm- 2 ] 


[keV] 


nt c Metal abundances' 1 Flux° 
[10 10 s _1 cm -3 ] [solar] [10~ 13 ergs cm -2 s] 


Lx f X 2 /d-o.f 
[10 34 crgs s- 1 ] 


13.9(10.7-17.2) 


6.1(3.1-25.8) 


1.7(1.3-2.7) 4.5(1.6-9.9) 8.2 


1.3 41.3/40 



Note. — Errors are at 90% confidence level. 

Note. — All parameters of the model (including normalization) are best-fitted with both ASCA and Chandra data points. 
a Hydrogcn column density. 

b Tcmperture of electrons in a thin thermal plasma, 
ionization parameter is an NEI. 

d Mctal abundances (helium fixed at cosmic). The elements included are C, N, O, Ne, Mg, Si, S, Ca, Fe, Ni. 
e Observcd flux in the 2.0-10.0 kcV band. 
'Absorption-corrected luminosity in the 2.0-10.0 keV band. 



Table 2 

Comparison of the parameters between GO. 570-0. 018 and SN 1987A 





Distance 


radius of an X-ray shell 


t 




[kpc] 


[PC] 


[years] 


G0.570-0.018 


8.5 


0.4 


80 


SN 1987A 


50 a 


0.12 a 
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a The parameters of SN 1987A are taken from Burrows et al. (2000) 
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